Systematics of radial and angular-momentum Regge trajectories of light non-strange qq-states* 
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We reanalyze the radial (n) and angular-momentum (J) Regge trajectories for all light-quark states with 
baryon number zero listed in the 201 1 edition of the Particle Data Tables. The parameters of the trajectories are 
obtained with linear regression, with weight of each resonance inversely proportional to its half-width squared, 
(r/2) . That way we are side-stepping possible channel-dependent and model-dependent extractions of the 
resonance parameters and are able to undertake an error analysis. The method complies to the fact that the pole 
position of the resonance is typically shifted from channel-dependent extractions by --^ r/2. This is also a feature 
of the large-Nc limit of QCD, where the masses change by r/2 when evolving from Nc = 3 to A'c = °o. Our value 
for the slope of the radial Regge trajectories is a = 1.35(4)GeV^. We discuss the fundamental issue whether 
the masses of the light-quark non-strange states fit into a universal pattern M„y = a{n + J)+b, as suggested by 
Afonin, and also predicted by some holographic models. Our joint linear-regression analysis in the {n,J,M ) 
Regge planes indicates, at a statistically significant level of 4.5 standard deviations, that the slopes of the radial 
Regge trajectories are larger from the angular-momentum slopes. Thus no strict universality of slopes occurs in 
the light non-strange meson spectra. 

PACS numbers: 14.40.-n, 12.38.-t, 12.39.Mk 
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I. INTRODUCTION 

The study of regularities in the hadronic spectrum has been 
a recurrent subject in the quark model 1 1 1, as it allows not only 
to check our current understanding of strong interactions, but 
also to predict possible missing states. In the case of light- 
quark mesons, which is the subject of the present study, the 
quark-hadron duality |2) implies QCD constraints based on 
the Operator Product Expansion of a two-point correlation 
function with some given mesonic quantum numbers (say J). 
In particular, 

f^j / (dM^^j /dn)^ const, (1) 

where M„j is the n-th mass of the meson and /„/ the cor- 
responding vacuum decay amplitude. More than a decade 
ago Anisovich, Anisovich, and Sarantsev |3| suggested that 
mesons could be grouped into radial Regge trajectories of the 
form 



M,? 



M^+nH^ 



(2) 



where Mq is the mass of the lowest-lying meson on each cor- 
responding trajectory and ji^ is the slope parameter Accord- 
ing to Ref. |i3J, the slope is approximately the same for all 
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the trajectories considered: jj.^ — 1.25(15) GeV^. The uncer- 
tainly was estimated based on the spread of the different re- 
sults for each meson family. In addition, some missing states 
predicted from Eq. (J2]) have indeed been confirmed and in- 
cluded in the latest edition of the Particle Data Group (PDG) 
tables |4|. Furthermore, Ref [3J also analyzed the venerable 
angular-momentum Regge trajectories (for a review see, 
e.g., i6j|), which motivated the original (rubber) string mod- 
els Q (for a review see, e.g., |8, 9|). Moreover, the large 
degeneracy ifTOl of the daughter Regge trajectories is capa- 
ble of producing the Hagedorn growth of the hadronic spec- 
tra ifTOinn (see Refs. 1 12, 13 1 for a recent reanalysis). 

In a remarkable paper Afonin llT4l (see also ifTSl ) analyzed 
jointly the radial and angular-momentum trajectories and ar- 
gued that they merge into a single pattern 

M^{n,J) = a{n+J)+c, (3) 

unveiling a kind of hydrogen-like accidental degeneracy, with 
a harmonic oscillator mass-squared spectrum. All these phe- 
nomenological findings provide some confidence on the string 
picture of hadrons, where the square of the mass is the fun- 
damental dynamical quantity. Together with the QCD short- 
distance constraint of Eq. (fill we may then infer that mesonic 
vacuum decay amplitudes tend to a constant in the upper part 
of the spectrum. 

Regardless of the success of the radial Regge trajectories, 
it is important to note that the resonance parameters, such as 
mass, width, or coupling constants, depend on the definitions 
and are sensitive to the background, i.e., to the particular pro- 
cess used to extract the resonance from the experimental data. 
This poses the relevant question of what the precise meaning 
of Eq. (|2| is, and, moreover, in what sense is QCD compati- 
ble with such an analysis. In the present work we reanalyze 



this problem, carrying out global linear regression fits with the 
uncertainty of the resonance position proportional to its width, 
r. Specifically, we use weights inversely proportional to the 
square of the resonance half-width. The approach is consis- 
tent with the fact that the pole position of the resonance is 
typically shifted from channel-dependent extractions by about 
r/2. Also, within the large-A^,. QCD [16', 17] (see e.g. |T8l 
for a review), where the strong coupling constant is assumed 
to scale as g ^ 1/ ^/Wc, the meson masses change by r/2 
when evolved from Nc = 3 to Nc = °°, as has been exploited 
intensely in Refs. |fT9lJ26<. 

We note that within the AdS/CFT proposal (for a review 
see, e.g., ||27l ) there have been attempts to formulate holo- 
graphic models (the so-called soft-wall models) with linear 
confinement |28| and, likewise, their light-cone relatives 1291 , 
complying to the ansatz of Eq. (J2]i. We recall that all these 
AdS/CFT inspired models are claimed to operate for large 
t'Hooft couplings, i.e., g ^ l/Nc- 

As we will elaborate in detail, our main finding, after 
considering the resonance width uncertainties, is to confirm 
the result of Ref. |3| with the updated data, as we find 
Ij} = 1.35(4) GeV^. On the other hand, our analysis in the 
{n,J,M^) Regge planes shows that at a statistically signifi- 
cant level of 4.5 standard deviations the slopes of the radial 
Regge trajectories are larger from the slopes of the angular- 
momentum trajectories. Therefore no strict universality of 
slopes occurs in the light non-strange meson sector. 

The plan of the paper is as follows. In Section III] we moti- 
vate our choice for the weight in the linear regression analysis. 

In Section [III] we discuss in detail, through the use of the 
present PDG tables, how the different states are grouped into 
the radial Regge trajectories. Whenever possible, we try to 
keep the successful choice of Ref. 1 3 1 taking into account the 
assumed uncertainties. In Section |III H| on we enlarge the 
choice of Ref. |3| to complete all the light unflavored states 
collected in the PDG. The update of the angular-momentum 
Regge trajectories is considered in Section IV In Section [V] 
we discuss, as originally suggested by Afonin, simultaneously 
the radial and angular-momentum trajectories. Finally, in Sec- 
tion |Vl] we summarize our results and draw our main conclu- 
sions. 

Throughout this work we use the up-to-date edition of the 
PDG tables Q. The symbol q stand for the light quarks, u or 
d. 



II. UNCERTAINTIES OF RESONANCE POSITIONS 

As already mentioned in the Introduction, in order to prop- 
erly size the meaning of the radial Regge trajectories for res- 
onant states it is important to review the well known fea- 
tures of the quantum mechanical decay process relevant to 
our discussion. The rigorous quantum-mechanical definition 
of a resonance with given quantum numbers corresponds to a 
pole in the second Riemann sheet in the (analytically contin- 
ued) partial-wave amplitude of the considered scattering chan- 
nel ll30l . This definition becomes independent on the back- 
ground, whereas the corresponding residue provides the am- 
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FIG. 1: The ratio of width to mass for the light-quark meson states. 
The surface of each point is proportional to the (27 + 1 ) spin degen- 
eracy, while the intensity is proportional to the isospin degeneracy 
{21+ 1). The band correspond to the average ± standard deviation 
bounds, r/M = 0.12(8). 



plitude to produce that resonance in the given process. 

However, although quoting the complex pole and the com- 
plex residue would be superior and highly desirable, for prac- 
tical reasons this is not what one typically finds in the PDG 
tables |4|, with very few exceptions. As a matter of fact, sev- 
eral definitions besides the pole in the second Riemann sheet 
are employed, such as a pole in the /T-matrix, the Breit- Wigner 
resonance, the location of a maximum in the speed plot, time 
delay, etc. (see, e.g., 113111321 ). 

A resonance may be interpreted as a superposition of states 
with a given mass distribution on the real axis, approximately 
spanning the M ± r/2 interval. Of course, the shape of the 
distribution depends on the particular process in which the 
resonance is produced, and thus on the background. Clearly, 
while all the definitions converge for narrow resonances, even 
for broad states we expect the masses obtained from various 
methods to be compatible within their corresponding M±r/2 
intervals. As stated above, the values listed by the PDG for a 
given resonance correspond to different choices of the def- 
inition and/or production processes, but mostly the results 
are compatible within the estimated width differences. This 
clearly provides an upper bound on the uncertainty of the res- 
onance position for different resonance parameter definitions. 
For shortness, we refer to this mass uncertainty estimate of the 
resonance mass as the half-width rule '. 

Quite remarkably, there is a QCD scenario where the half- 
width rule estimate becomes parametrically small for all the 
resonances in the mesonic spectrum. In the large-A^^ limit 
of QCD lfT6l |T7 | (see e.g. |18| for a review using effec- 
tive Lagrangians) mesons become stable, i.e., their masses are 
M = ^iN^), while their widths are suppressed, F = ff{l/Nc), 
such that the ratio Y /M = ff{\/Nc). This expectation of the 
large-A^c limit seems to be fulfilled very well in the real Nc = 3 



' Of course, the width itself has an uncertainty which may eventually enlarge 
the global indetermination in the resonance mass. 



world, since one finds for the light-quark mesons an average 
value r/M — 0.12(8) (to be compared with a rule-of-thumb 
l/Nc = 0.33 for A^^ ^ 3) |l33l. This feature is visuaHzed in 
Fig. [T] Of course, there are exceptions to this average ra- 
tio, but they are scarce within the given confidence interval. 
In fact, only just one state {a = /o(600)) goes over the 1/3- 
value f33l. 

A complementary way of connecting parametrically the 
mass shift and the decay width is as follows. One starts with 
the leading-A^t resonance Lagrangian |18|, recalling that the 



three- and higher n-mesonic interactions are ^(A^, 



l-n/2. 



m 



[TtI . Thus, the mass shift is computed as a loop integral via 
the self-energy whose imaginary part corresponds to the de- 
cay width of the particles inside the loop according to the 
Cutkosky rules. This argument makes it clear that the l/Nc 
scaling of the mass-shift and the width are exactly the same 
and bound by &{\/Nc), although the numerical values of the 
two quantities may not coincide exactly. The point of this dis- 
cussion is that if we take the leading- l/A^^ resonance mass, 
its systematic uncertainty is parametrically indistinguishable 
from the decay width, since they are the real and imagi- 
nary parts of the self-energy, respectively. As pointed out in 
Ref. 1341 . the role of the mass-shift is crucial when determin- 
ing the properties of two-point correlator functions. 

Within this framework, the half-width rule has been used 
recently ll33l [35l for the case of the scalar and pseudoscalar 
mesons with rather interesting results regarding the identifica- 
tion of glueball states and chiral symmetry doublets. Here we 
extend these ideas to the rest of the light-quark meson spec- 
trum. Specifically, to incorporate the half-width rule in prac- 
tice, we take ^ 



x' = 



Mi-M: 



H.exp 



r„M„ 



(4) 



for the linear regression fit, where the radial Regge formula, 
Eq. p), is used as the model. Note that in doing so, we are 
just saying that Eq. S is fulfilled within the uncertainty M^ = 
jl^n +Mq ±r„M„. Moreover, we will stay within the linear 
ansatz as the half- width rule yields insensitivity to small non- 
linearities as analyzed e.g. in Ref. f36l for n-trajectories or in 
Ref. Ii37il for J-trajectories. 



III. RADIAL REGGE TRAJECTORIES 

The construction of a meson Regge trajectory requires a 
choice on the possible meson assignments. The analysis of the 
radial Regge trajectories we are candying out consists of two 



which does not alter 



■ There is an alternative fit with X^ = Y.n\ "r z""'' 
much in the results. Actually, both x^ functions are particular examples of 
the more general maximum-likelihood method, where the resonance pro- 
duction profile is assumed to be Gaussian. For a discussion on other pro- 
files, in particular for the ubiquitous Breit-Wigner shape, see Appendix [a] 
for details. 




FIG. 2: (color online) The {n,M ) plots for the states ai (1++) (lower 
solid and dashed lines) and ai(3++) (upper solid and dot-dashed 
lines). Error bars correspond to take AM = iPM. 



stages: The first one reanalyzes the results of Ref. ||3l with 
the inclusion of more states from the updated PDG tables r4l, 
while from Sec. lIIIHl on we deal with meson families not con- 
sidered in Ref. O. To facilitate the comparison, we follow as 
close as possible the presentation of Ref. |3 1. 

We motivate our selections with rather detailed discus- 
sions. The reader interested in the results only may jump to 
Sec. inrN] 

In all our M^ -plots we take, in line with the half-width rule, 
the error to be given by AM^ = iPM. 



A. fli (11++) and 03(13++) 

Compared to Ref. Q, we consider four different trajec- 
tories: two for the ai(ll++) states (lower solid and dashed 
lines in Fig. Kb and two for the 03(13++) states (upper solid 
and dot-dashed lines in Fig. |2|i. The first trajectory for the 
ai(ll++) states contains fli(1260), fli(1640), and fli(2095). 
The ai (2340) state (now called oi (2270)), assumed in Ref. ^ 
to belong to this trajectory, is now used in the daughter trajec- 
tory for the ai(ll++), together with a new state not consid- 
ered in Ref. ||3|, the ai(1930). The linear fit to the first trajec- 
tory for the fli(ll++) states yields /i^ == 1.36(49)GeV^ with 
j;i^^/DOF ==0.12 (lower-solid line in Fig. pj). For the second 
trajectory n^ = 1.43(73) GeV^ (dashed Une on Fig.|2|. This 
trajectory has only two states and it will not be considered for 
the final compilation. 

In the case the a3(13++) trajectories, the first one contains 
the new a3(1875) state together with the a-i(2275), yielding 
/i^ = 1 .5(1 . 1) GeV^ (upper solid line in Fig.pb and the second 
trajectory contains only the 03(2030) state (dot-dashed line 
parallel to the 03(13++) trajectory). 
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FIG. 3: (color online) The {n,M'^) plots for the r}(00"+) (lower solid 
and dashed lines) and T)2(02 '") (upper solid and dashed lines) tra- 
jectories. Error bars correspond to take AM = ±TM. 



FIG. 4: (color online) The {n.M^) plots for the states pi(ll ) 
(solid and dashed lines, respectively), and the p3(13 ) states (dot- 
dashed line). Error bars correspond to take AM = iFM. 



B. T](00-+)andT]2(02- 



C. pi(ir-)andp3(13- 



Figure p] shows the 77(00 ^) and Tj2(02^+) states where, 
due to two independent flavor components qq and ss, both 
yield two trajectories. The 77(00 ^)qii contains five states: 
77(548), 77(1295), 77(1760), 77(2100), and 77(2320) (lower 
solid line in Fig. Isll, where the first state, 77(548), is not used 
in the linear fit. In Ref. [3] the state 77(2100) was predicted, 
while nothing was said about the 77 (2320). Both states, now 
listed by the PDG, are incorporated in our study. The fit yields 
li^ = 1.33(1 1) GeV2 with x^ /ViOV = 0.26. 

The 77(00"+)ij trajectory with four states: 77(958), 
77(1475), 77(2010), and 77(2225), yields ^^ = 1.36(14) GeV^ 
with xVDOF = 0.44. The 77 (2010) was predicted in Ref. |3 1 
under the name 77 (1900) and now is listed in the PDG Tables. 

In Ref. 131 only one state with mass near 1440 MeV was 
considered. Now it is well established that in this energy re- 
gion there are two different 77 states, the 77(1405) and the 
77(1475). The first one, however, is not unambiguously lo- 
cated and it is considered to be a glueball (see the mini-review 
about this state on the PDG Tables), therefore we exclude it 
from our fitting procedure. The second state is included in the 
ss trajectory. 

The 77 (02^+),,,^ trajectory yields jx^ = 1.32(32) GeV^, 
with x^/DOF = 0.22. This trajectory contains 772(1645), 
772(2030), and 772(2250). The 77(02"+).vs trajectory, which 
contains only one 772(1870) state, is drawn parallel to the non- 
strange case. 



The two trajectories for pi(ll ) are depicted in Fig. 
The first one contains p(770), p(1450), p(1900), and 
p(2150). As explained by the PDG, it is not clear what val- 
ues for the mass and width one should use for p(1900). We 
choose M = 1 .870(30) GeV and F = 0. 150(20) GeV. The Hn- 
ear fit (solid line in Fig. [ijl yields /z^ = 1.43(13) GeV^ with 
xVDOF = 0.09. 

The second trajectory contains p(1700), p(2000), and 
p(2270). These last two states were predicted by Ref. lO 
and now are listed in the PDG compilation. The p(1700) 
and p(2000) states, however, are controversial and need 
confirmation. The corresponding slope trajectory is /i^ = 
1.08(47) GeV^ with a jVDOF = 0.004 although it is drawn 
in Fig. |4] as parallel to the qq trajectory due to the lack of 
confirmation of these states. There is a new state in the 
PDG Tables called p(I570), which also needs further con- 
firmation because it might reflect a threshold effect or an OZI 
suppressed decay mode of the p(1700) (see the mini -review 
about this issue on PDG). We do not include it on our analysis, 
either. 

In addition to Ref. f3l, we have also considered the 
P3(13"") states, which include p3(1690), P3(1990), and 
P3(2250). The slope for this trajectory is ^2^ 2 29(32) GeV^ 
with ;t;VDOF = 0.05. Neither P3(1690) nor P3(2250) are 
well established resonances and we just quote them for com- 
pleteness. We do not use the slope prediction for this trajec- 
tory in our later average result. 
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FIG. 5: (color online) The (n.M^) plots for the ;r(10"+) (lower solid 
line) and 7^2(12 '") (upper solid and dashed lines) trajectories. Error 
bars correspond to take AM = iPM. 
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FIG. 6: (color online) The (n^M^) plots for the ao(10++) (lower 
solid line), a2{yi'^^) qq and ss (upper solid and dashed lines, re- 
spectively), and 04(14++) (dot-dashed line) trajectories. Error bars 
correspond to take AM^ = iPM. 



D. K and K2 

For the ;r(10 ^) trajectory (lower solid line in Fig.pji, com- 
posed of ;r(1300), ;?r(1800), ;?r(2070), and ;r(2360)7the fit 
produces li^ = 1.27(27)GeV2 with ;t;2/DOF = 0.16, where 
the stable 7i:(140) state is not used in the fit ^ . We update the 
trajectory including the ;r(2070) and 7i:(2360) states originally 
predicted in Ref . 1 3 1 . 

The 7^2(12^+) states produce two trajectories. The first one 
includes 712(1670) and two new states predicted in Ref. lIS): 
7:2(2005) and 7r2(2285). The fit yields ^^ ^ 1.21(36) GeV^ 
with ;if^/DOF ~ 0.02 (solid-upper line in Fig. pb. The sec- 
ond trajectory has two new states, not predicted in Ref. |l3l: 
7r2(1880) and 7r2(2100). The heaviest 7r2, with the mass of 
2.090(29) GeV and the width of 0.625(50) GeV, has still to 
be confirmed. Conversely, if we use the fitted daughter trajec- 
tory with this state omitted, we predict its mass to be around 
2.19(13) GeV. 



E. ao(10++), 02(12++), and 04(14++) 

In Ref. IS the experimental information in the flo(10++), 
(22(12++), and 04(14++) sector was scarce and could not fix 
the /i^ slope uniquely. Therefore, two different slopes, jti^ = 
1.38 GeV^ and jU^ = 1.1 GeV^ where deduced depending on 
the states included, and in fact jX^ = l.l GeV^ predicted a yet 



unobserved new state flo(1800). Hence /i^ ~ 1.38 GeV^ is 
favored currently and we accept the classification of Ref. |3|, 
with the ao(10^^), 02(12++), and 04(14++) trajectories and 
the 02(12^^) split into two daughters. The flo(10^^) con- 
tains flo(980), flo(1450)'*, and ao(2020). The prediction for 
the slope is pL^ = 1.42(26) GeV^ with X^ /XiO¥ = 0.48. The 
flo(2260) predicted in Ref. pl has not been seen yet and in our 
present description should be located around 2.29(12) GeV. 

Two trajectories for 02(12^^) are presented. The lower tra- 
jectory contains 02(1320), 02(1700), and 02(2175)^, giving 
ju2 = 1.39(26) GeV2 with a x"- /XiO¥ = 0.24. The previously 
predicted 02 (2400) Ref. 1 3 1 has not been seen yet. It is also 
predicted within our trajectory to have M = 2.42(17) GeV. 

The upper trajectory contains 02 (2030) and 02 (2255). The 
02(2030) is an average of different experimental determina- 
tions (under two different names) from the PDG compila- 
tion. In the 1999 PDG edition a state called 02 (1990) was 
introduced, while in 2001 this state was updated to become 
02(2030) by Ref. |38|, but not modified in the PDG review. 
Since then 02 (2030) appears under two different entries in the 
PDG compilation, hence one of them is redundant. The mass 
and width for 02 (2030) are not averaged by the PDG, where 
just the three different measurements are presented. We av- 
erage them with the result 2021(14) MeV for the mass and 
220(23) MeV for the width. This trajectory would produce 
^2 = l.0(7)GeV2. 

Another problem to face for the upper 02 trajectory is the 



^ One expects a strong non-linearity for the Goldstone bosons, see, e.g., 
Ref. PI. 



* Called before flo( 1520). 

'' Named before 02(1660) and 02(2100), respectively. 



presence of two very close resonances, 02(1950) with the 
mass of 1950(50) MeV and the width of 187(50) MeV, and 
02(2030). It is argued in Ref. |38| that it is necessary to ob- 
tain a better fit to the data. As a mater of fact, due to the large 
errors of the mass position and widths, these states might eas- 
ily be a single state. For the presented reasons, we do not use 
this trajectory for our average slope value. 

Finally, the 04(14++) trajectory contains O4(2040) and 
04(2255) and the slope turns out to be ji^ = 1.0(8) GeV^. We 
draw, however, a parallel line to the main (solid) trajectory in 
Fig.|6] 



values for both masses and widths, it might turn out that 
/2(1810) and /2(2220) are the very same /2(1910) and 
/2(2240) states. 

• The upper trajectory is described with ji^ = 
1.43(83) GeV^ and contains two states, /2(2010) 
and /2 (2340). 



G. /o(00++) 



F. /2(02++) 

In Ref l3\ it was not possible to discriminate the slopes 
^2 = 1.1 GeV2 or jx^ = 1.38 GeV^ for the /2(02++) sector 
with the 12 available states in year 2000. The currently listed 
18 states favor the second slope. 

Figure 2(c) of Ref. ID shows a quadruplet of trajectories 
wtih two flavor components, qq and ss. With the inclusion of 
the additional 6 new states, we find that an overall satisfactory 
update of Ref. |3] is given by the scheme presented in Fig. [7] 
(left panel) requiring some reshuffling which we describe be- 
low. We name these trajectories, /2 ,/2 ^ /^ and /|'. 

• The /j trajectory (lower solid line in Fig. T\. It contains 
/2(1270), /2(1750), and /2(2150), one less state than 
Ref. f3l which also included /2(2400). In our case, the 
slope for this trajectory is /i^ — 1.50(19) GeV^ with 
xVDOF = 0.06. 

• The /I trajectory (dashed line in Fig. rTb. It contains 
/2(1430), which is still to be determined (we take M — 
1468(60) MeV and r= 100(100) MeV),/2(1910) with 
M = 1927(32) MeV and F = 154(73) MeV, and, finally, 
/2(2240). These states yield ji^ ^ 1.48(23) GeV^ with 
xVDOF = 0.09. 

• The /j trajectory (the upper solid line in Fig. iTl) As 
in Ref. 13|, it is composed of /2(1525), /2(1950), and 
/2(2295), giving jx^ = 1.47(25) GeV^ with jVDOF == 
0.00001. 

• The /j' trajectory (the dot-dashed line in Fig. rTb. It con- 
tains /2(1565) (which needs confirmation), /2(2000), 
and /2(2300). The slope for this trajectory is jj.^ ~ 
1.42(20) GeV2 with jVDOF = 0.05. The states 
considered here involve some reshuffling compared to 
Ref. 0. 

We now turn to the new trajectories, i.e., not given in 
Ref. O, which are separately plotted in Fig. IT] (right panel). 

• The two trajectories including /2(1640) and /2(2150) 
as weU as /2( 18 10) and /2(2220) (both of them need 
confirmation), might actually be intertwined or describe 
an overcomplete set of states. The first one returns /i" = 
1.99(36) GeV2 and the second [i^ = 1.69(36) GeV^. 
Considering the lack of confirmation and the particular 



Two trajectories for /o(00++) are displayed in Fig.mland as 
claimed in Ref. 1 3 1 they are doubled due to two flavor compo- 
nents, qq and ss. Without considering the /o(600) (see how- 
ever Ref. 1 35 1 and Appendix [A)i, also called the a meson, 
the lower trajectory contains four states: /o(980), /o(1500), 
/o(2020), and /o(2200) (solid fine in Fig.js]). The last state 
was actually predicted in Ref. f3| and later confirmed exper- 
imentally. The trajectory yields \i} — 1.31(12) GeV^ with 

xVdof = o.ii. 

The second trajectory (dashed line in Fig. [8]) has also four 
states, /o(1370), /o(1710), /o(2100), and /o(2330), where an 
average of the experimental determinations is considered for 
this latter state. It yields ^^ ^ 1.24(18) GeV^ with;t;VDOF = 
0.12. 



H. Q)(01--)andco3(03"") 

After reanalyzing the radial Regge trajectories taken into 
account in Ref |3 1, we now analyze using the same methodol- 
ogy the remaining meson families included in the latest PDG 
review [4|. 

Two trajectories for the a)(01 ) states are shown in Fig.^ 
The ordering of the states on the different trajectories for 
the o-family follows very closely the classification of the p- 
family states, Fig.ffl The main o-trajectory, representing the 
qq states and drawn as a solid line in Fig. l9j contains four 
states: CO (782), CO (1420), CO (1960), and co (2205). The slope 
from the fit is \i^ = 1.50(12) GeV^ with ;t;VDOF = 0.32. 

The daughter trajectory (dashed line in Fig.l9]l contains two 
states: CO (1650) and CO (2290). In the PDG recollection, two 
(i) states with very similar masses are listed in the 2.3 GeV 
region: co(2290) and co(2330). Looking at the error deter- 
mination of the parameters of these states it is not clear to us 
that the two states are indeed different. We gather both ex- 
perimental results in a single entry co (2290), which has mass 
2315(45) MeV and width 325(185) MeV. Comparing this tra- 
jectory with the corresponding one from the p -family, we no- 
tice a missing (i) state with the mass near 2000 MeV, indicating 
the state CO (2290) to be the third on its trajectory. With only 
two states, the corresponding slope is [i^ = 1.27(47) GeV^. 

The third trajectory in Fig. ^ describing the C03(03 ) 
states, contains C03(1670), C03(1945), and C03(2255). It yields 
H^ = 1.16(26) GeV2 with ;t;VDOF = 0.37. 
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FIG. 7: (color online) The {n,M ) plots for the seven /2(02++) trajectories comprising the 4 originally described in Ref. |3l (left panel) and 
the new ones discussed in the main text (right panel). Error bars correspond to take AM^ = iPM. 
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FIG. 8: (color online) The (n.M^) plot for the /o(00++) qq (solid 
line) and .si^ (dashed line) trajectories. Error bars correspond to take 

am2 = ±rM. 
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FIG. 9: (color online) The {n,M^) plot for the co(01 ) (circles, 
solid line and triangles, dashed line) trajectories, as well as for the 
0)3(03^^) (squares, dot-dashed line) trajectory. Error bars corre- 
spond to take AM^ = iFM. 



I. /ii(01+-)and/i3(03+-) 

The h sector contains two trajectories corresponding to the 
/ii(01+_) and /!3(03+^) states shown as circles and squares 
in Fig. 10 respectively. The /!i(01+^) case consists of four 
states: /!i(1170), /ii(1595), /!i(1965), and /ii(2215). The 
/ii (1380) state is excluded, since it still needs to be confirmed. 
The linear fit to this trajectory, shown as a solid line in Fig. 10 



gives ^2 = 1.20(25) GeV^ with jVDOF = 0.01. 

The /!3(03+^) includes only two states, /!3(2025) and 
hj,{2275) (supposed to be the second and third excitation 
states of that trajectory), thus the slope /z^ = 1.08(54) GeV^ 
is determined. In Fig. 10 the dashed line, drawn parallel to 
the solid line, represents this trajectory. 
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FIG. 10: (color online) The {n,M^) plot for the /ii(01+") (circles, 
solid line) and /i3(03^ ) (squares, dashed line) trajectories. Error 
bars correspond to take AM = iPM. 
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FIG. 12: (color online) The («,M^) plot for the /i(01++) (circles, 
solid line) and /3(03++) (squares, dashed line) trajectories. Error 
bars correspond to take AM = iFM. 
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FIG. 11: (color online) The (n.M^) plot for the foi(ll+") (circles, 
solid line) and ^3(13^ ) (squares, dashed line) trajectories. Error 
bars correspond to take AM = ±TM. 



J. /?i(01+-)and&3(03+-) 

Similarly to the previous subsection, the b sector contains 
two trajectories, coiTesponding to the /?! ( 1 1 +^ ) and ^3 ( 1 3+^ ) 
states, shown as circles and squares in Fig. [TT] respectively. 
The fei(01+") consists of three states: bi{\235), b\{\96Q), 
and ^71 (2240). The linear fit returns jU^ = 1.17(18) GeV^ with 
;t;VDOF = 0.00001 (soHd line in Fig. [IT 1. 



The ^73(03+ ) includes only two states, ^73(2030) and 



^73(2245), hence the slope is pi- = 0.93(75) GeV^. In Fig.[ll 
a dashed line, parallel to the solid line, represents this trajec 
tory. Since the resemblance between the h sector and the b 
sector is apparent, that suggests the existence of a still not de- 
termined b\ state with a mass of the order of 1600 MeV. 



K. /i(01++)and/3(03++) 



The situation with the /i and f^ states is equivalent to 
the b and h case, we thus have two different trajectories 
coiTesponding to the different angular-momentum, /i(01++) 

re- 



12 



and /3(03++), shown as circles and squares in Fig. 
spectively. The /i(01++) trajectory consists of three states: 
/i(1285), /i(1960), and /i(2240). The fit for this trajectory 
returns iT- = 1.19(15) GeV^ with X^ /T)0¥ = 0.13 and it is 
shown as a solid line in Fig. [12] 



The /3(03++) includes 
/3(2300). The slope is ^2 ^ 



only two states, f^ (2050) and 
= 1.27(64) GeV2 



12 



the 



In Fig. 

dashed line, drawn parallel to the solid line, displays this tra- 
jectory. The location of these states follows closely the case of 
/!3(03+^) and ^73(13+^), hence it starts at the radial quantum 
number n = 2. From Fig. 12 it is not obvious how to allo- 



cate the /i (1420) state, since the departure from the expected 
value seems much larger than expected from the half-width 
rule, therefore we exclude it from the fit. This choice resem- 
bles the h\ case. 
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FIG. 13: (color online) The {n,M^) plot for the 0(01 ) (circles, 
solid line) and 03(03^^) (squares, dashed line) trajectories. Error 
bars correspond to take AM = iPM. 
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FIG. 14: (color online) The {n,M'^) plot for the ;ri(ll +) (circles, 
solid line) trajectories. Error bars correspond to take AM^ — iFM. 



L. il){0i — )andit)i{03—) 

The0 sector has three states with 7= 1 (0(1020), 0(1680), 
and (j) (2170)), and one with 7 = 3 (03 (1850)), therefore the al- 
location of states becomes less unique. However, if the states 
with 7 = 1 are placed on a radial linear trajectory, a well de- 
termined slope n^ = 1.84(6) GeV^ with a x^ = 0.06 is ob- 
tained. We note that such slope is much larger than any of 
the other slopes found so far. We could also concede the 
states 0(1680), and 0(2170) to ben = 2,3 respectively, which 
would produce /i^ = 1 . 19(4) GeV^ instead, although with too 
large ;if^/DOF = 6.4. Due to this ambiguity, we will not con- 
sider this family for the final summary results 

In Fig. 13 two trajectories are shown as parallel lines, with 
a solid line representing the 7=1 states, and the dashed line 
going across the the 7 = 3 state. Clearly, this somewhat dis- 
turbing picture should profit from both theoretical or experi- 
mental insight. 



M. Ki{ll-+) 

Finally, the last sector we analyze corresponds to the TZi 
states, composed by ;ri(1400), ;ri(1600), and ;ri(2015). Two 
different measurements are found for this last state and we 
average them to have mass M = 2.013(25) GeV and width 
F = 0.287(53) GeV. 

With these three states on a linear trajectory the slope ob- 
tained is n^ = 1.09(36) GeV^ with a ;(;2 = 0.11. The results 
are depicted in Fig. [14] 



N. Summary of the radial-trajectory fits 

We summarize this Section by collecting the fits for all the 
radial trajectories studied. The jx^ parameter ranges from 
1.09(36) GeV2 (the ;ri(ll-+) trajectory) to 1.50(19) GeV^ 
(corresponding to the /I (02++) trajectory). The weighted av- 
erage yields ^ 



jLi^ = 1.35(4) GeV^ 



(5) 



This result agrees within the uncertainties with the estimate 
of Ref. (3], 11^ = 1.25(15) GeV^, where the uncertainty was 
given by the spread of the mean values determined from a 
fit to the PDG masses with equal weights. When we carry 
out the same procedure for the updated and new trajectories, 
we get jx^ = 1.32(12) GeV^. This seems to provide a quite 
robust estimate of a common radial Regge trajectories slopes. 
A graphic overview of the estimated slopes is presented in 
Fig. [T6] (upper part). 

We have also considered the possibility of a linear n- 
dependence of the masses, since it was a popular outcome 
of holographic models in the hard-wall scheme (see e.g. 
Ref. 1221 and rererences therein). With the same conditions as 
analyzed above, i.e., assuming the validity of the half-width 
rule our analysis is not compatible with such radial spectrum; 
typically we obtain ;i;^/DOF ~' 10 or larger. 



^ We use the customary definition for the weighted average A = 
Y!i=i i^i^i/EiLi >''/. with Wj = 1/of . The eiTors are the mean-squai'ed de- 
viation, y^A^ — (A)^. The mean average corresponds to vy, = 1 . 
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X 


m|(0) [GeV2] 


fe [GeV2] 


a^(0) [GeV-2] 


X^/DOF 


n 


0.30(0) 


1.22(10) 


0.82(7) 


0.37 


p 


0.60(11) 


1.19(10) 


0.84(7) 


0.15 


n 


0.018(0) 


1.29(11) 


0.78(7) 


0.26 


ai 


-0.45(43) 


1.09(18) 


0.92(15) 


0.03 


ao 


0.96(7) 


1.02(12) 


0.98(12) 


0.001 


h 


-0.66(48) 


1.15(16) 


0.87(12) 


0.12 


h 


2.26(16) 


1.20(17) 


0.83(12) 


0.01 


f. 


0.96(7) 


1.13(12) 


0.88(9) 


0.07 


CO 


-0.48(11) 


1.09(11) 


0.92(9) 


0.02 



TABLE I: The {J,M ) trajectories for leading and daughter trajec- 
tories. For an easy of comparison with Ref. |3|, we also show the 
corresponding 0!^(0) ~ l/jS for each trajectory. 



IV. {J,M^) TRAJECTORffiS 

In this section, taking into account all the states so far con- 
sidered and adding those with larger 7 = 4,5,6 from the PDG 
Tables |4|, we complement the results of the {n,M^) analysis 
with the study in the {J,M^) plane, i.e, the standard Chew- 
Frautschi plots Q . 

One may parametrize the trajectories as 



with (Xx{0) and 0:^(0) constant parameters. 
Equivalently, we consider 

M|(y)=M|(0)+j3x/. 



(6) 



(7) 



The j3 parameter is related to a as j3 ^ l/o:^(0). 

The results for the leading trajectories are shown in Tablell] 
Generally, we do not attempt to determine the slope when the 
trajectory is made of less than three states. An overview of the 



estimated slopes can be seen in Fig. 16 (lower part). Never- 
theless, to provide a broader perspective, we show in Fig. 15 
also trajectories with just two states (dashed lines). 

The weighted average result for the angular trajectories 
yields 



j3 = 1.16(4) GeV^ 



(8) 



If one considers, instead, the spread of central values as in 
Ref |3l, the updated result for the states, one obtains j3 = 
1.15(8) GeV^, in agreement with |3|. 

When comparing the results of Eq. (J5]l and Eq. ([8]), we note 
that the radial and angular-momentum slopes are different at 
the level of 3.4 standard deviations. Thus, there is indication 
that the radial slopes are larger than the angular-momentum 
slopes at a significant statistical level. We come back to this 
important issue in the next section. 

It should be noted that in addition to the states of Ref. |f3l, 
we have also included the O), hy, bi, /i, (j), and tti sectors on 
our analysis, as well as three new 7-trajectories, with the ao, 
/o, and (0 states. 



V. JOINT {n,J,M^) FITS 



A summary of the estimated radial and angular slopes is 



given in Fig. 16 We consider a weighted average estimate for 
a common (universal) slope trajectory, including the studied 
trajectories both in the {n,M^) and the {J,M^) planes. This 
yields /z^ = j3 = 1.26(3) GeV^. When the spread of central 
values for all the trajectories is used, we obtain /i^ = j3 = 
1.27(14) GeV^. The trend to produce a number closer to the 
radial slope reflects a larger sample. As a matter of fact, when 
we weight both n and J trajectories equally, we get jti^ = ^ = 
1.26(3) GeV2. 

In the previous sections we have caiTied out the analysis for 
the {n,M^) and {J,M^) planes, with the indication that /X^ > j3 
at a significant statistical level (3.4 standard deviations). This 
strongly suggests a careful reconsideration of the findings of 
Ref. lfT4ll . where a common fit of Eq. (pji for the («,M^) and 
(y,M^) planes was proposed. Technically, the joint analysis 
presented in this Section is different form the separate analy- 
ses of Sect. Ill and IV in the following important detail. In 



9 /Z 

the separate fits the constants Mq and M q in the formulas 



^0 



M'q were treated as un- 



M^{n)^fi^n+Ml^ andM2(«) : 
related parameters, even in the same family of states. On the 
contrary, formula (Isl), with a common parameter c for a given 
family, relates the "offset" constants Mq and M q, providing a 
constraint to the statistical analysis. 

Therefore, to look closer at the issue of universality of the 
radial and angular-momentum slopes, we analyze each sector 
independently with two different plane fit functions: the non- 
universal formula 



M =an + bJ + c 



on one hand, and the universal formula 



M^ ^a{n + J)+c 



(9) 



(10) 



17 for the p- and 



on the other hand. With the c parameter fixed for the whole 
family, the states with different values of n and J belong to 
a Regge plane. This is illustrated in Fig. 
Tj -families as an example. 

Our numerical results are collected in Table [ill Several me- 
son sectors can be placed on two almost parallel planes. The 
nomenclature used is as follows: after the name of each fam- 
ily, the subindex quoted refers to the states with the particu- 
lar angular quantum number used on the plane. For example, 
Pi35 means the set of all the p states with angular-momentum 
J = 1,3,5, and with all the possible radial quantum numbers. 

The numbers in Table III] show a few interesting features. 
The most important one is that for each family of states the fit 
with Eq. (J9]l is preferred over the fit with Eq. ( 10 1 (judging by 
the different ;if^/DOF values). Moreover, we generically find 
a > b. The offset parameter c also seems to be stable through 
all the planes, although less stable than the radial and angular 
slope parameters. 

The fact that the ratio of the radial to angular slope, denoted 
as R, grows with the quark mass may be a generic and phys- 
ically relevant feature. Note that for the heavy quarkonia the 
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FIG. 15: (color online) The (J,M~) plots for all states considered in Tableyj Error bars coiTespond to take AM = iPM. 
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FIG. 16: (color online) (n.M ) and {J.M ) slope results for the con- 
sidered trajectories. The horizontal dashed line separates the radial 
slopes (circles) from the angular slopes (squares). Individual errors 
are estimated from the x fits to the corresponding trajectories de- 
scribed in the main text. The bands correspond to the weighted av- 
erages of the radial (upper band) and the angular-momentum (lower 
band). 



joint fit is compatible with the formula fvfi — a{2n + J)+c 
|p9l. Thus, it may be that R is close to unity for light mesons, 
equal 2 for heavy mesons, and assumes an intermediate value 
for hidden-strangeness states. 

From Tablelnlit is worth stressing how close the states from 
the p and the co families are. For 002 plane, we have assigned 
to 02(2030) and 02(2255) radial quantum numbers n = 2 and 
3, respectively. Otherwise (with n — 1 and 2), the angular 
slope would be larger than the radial slope. The /113, fen, and 
/i3 have systematically smaller a and b parameters from the 
remaining families, although very similar among themselves. 
As commented already, this fact may be caused by the lack of 
states in these sectors. States with higher angular-momentum 
(when discovered) would lead to better and more reliable de- 
termination of the Regge plane parameters. Similar comments 
apply to the ^13 plane. 

Considering only the planes with no large hidden 
strangeness content and with six or more states (excluding 
then the /113, bn, /13, and 0i3 planes)^ we obtain our global 
fit with the result 



Therefore, the a — jx^ parameter reads a — 1.38(4) GeV^ for 
the global Regge-plane fit, compatible with Eq. (Bl. The b — fi 
parameter reads b = 1.12(4) GeV^, also close to the value 
of Eq. dSll. We note that a > b at the level of 4.5 standard 
deviations. In this estimate we take the geometric average of 
the individual errors (equal 0.04) for the standard deviation of 
the difference a — b. Therefore the joint analysis points at a 
lack of universality of the Regge slopes. 

The right part of Table III] shows the result of the fit, where 
universality is imposed. This fit cannot be statistically rejected 
based on the values of ;if^, however, it is somewhat worse than 
without the universaUty constraint. 



VI. CONCLUSIONS 

In this paper we have reanalyzed, with the help of the up-to- 
date PDG tables [4], the linear radial and angular-momentum 
Regge trajectories considered in Ref. t3J, including in the fits 
the width of each state as an estimate of the error of the res- 
onance mass (the half-width rule). As we have explained this 
is a reasonable way to smooth out resonance profile informa- 
tion, which makes the very definition of the resonance mass 
ambiguous. Moreover, this choice allows to undertake an er- 
ror analysis, not carried out in Ref |3|. Furthermore, we have 
argued that such a procedure fully complies to the large-A^^ 
viewpoint and actually suggests an interesting interpretation: 
the Regge-fitted masses are considered to be the leading-A^^ 
contribution to the mass of the resonance. This incorporates 
a desirable flexibility as to what should the Regge fit be com- 
pared to. The squared mass of each meson is then represented 
as M^, = M^ ± FM, where M is its mass and F is its width. 

Generally, we reproduce the results of Ref. [3J when no 
uncertainties are included. This only reflects the robustness 
of the main features of the PDG compilation along the last 
10 years, although some numerical values of the masses have 
changed and, furthermore, some new states, partly predicted 
by the pioneering radial Regge analysis of Ref |3 1, have been 
added. From our results it follows that there is no need to 
consider further new states to get an acceptable Regge de- 
scription. This is consistent with an assertion of a complete 
mesonic spectrum up to the highest energies considered in our 
work. 

We have also addressed the issue of the universality of ra- 
dial and angular-momentum slopes within the errors deduced 
from the linear regression analysis with weights provided with 
the half- width rule. Our joint analysis in the {n,J,M^) Regge 
planes indicates, at a statistically significant level of 4.5 stan- 
dard deviations, that the radial slope is larger from the angular- 
momentum slope. Thus no strict universality of slopes occurs 
in the light non-strange meson spectra. 



M^ 



1.38(4)« + 1.12(4)7- 1.25(4). 



(11) 



Appendix A: Dependence on resonance profile 



^ The /02 is not consider either due to the arbitrariness on the selection of its 
components. 



In this appendix we show the independence of our results on 
the shape of the resonance profiles, and hence support our x^- 
statistical treatment. The ;if^-nature of the fit relies implicitly 



13 




FIG. 17: (color online) The (n,J,M^) Regge planes for the r\ family and the p family. 
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'702,, 


1.36(14) 


1.27(22) 


-0.50(28) 


0.44 


1.34(13) 


-0.47(27) 


0.35 


Pl35 


1.36(12) 


1.12(9) 


-1.87(23) 


0.40 


1.21(7) 


-1.79(22) 


0.73 


T^2A 


1.47(10) 


1.27(10) 


-1.45(10) 


0.29 


1.36(6) 


-1.34(6) 


0.50 


"246 


1.35(25) 


1.06(16) 


-1.75(48) 


0.14 


1.15(13) 


-1.70(48) 


0.33 


"02 


1.35(24) 


0.78(24) 


-0.39(27) 


0.53 


1.06(9) 


-0.09(13) 


0.90 


hlAd 


1.38(13) 


0.64(8) 


0.04(33) 


0.85 


0.76(8) 


1.06(29) 


5.03 


hi 


1.34(11) 


0.69(6) 


-0.38(15) 


0.14 


0.84(6) 


0.13(11) 


5.66 


«135 


1.42(11) 


0.98(8) 


-1.78(11) 


0.63 


1.16(5) 


-1.70(10) 


1.77 


hn 


1.17(23) 


0.75(19) 


-0.53(62) 


0.02 


0.93(13) 


-0.37(61) 


0.46 


bi3 


1.15(17) 


0.72(15) 


-0.35(32) 


0.05 


0.91(9) 


-0.29(32) 


0.88 


M 


1.19(15) 


0.70(19) 


-0.24(17) 


0.07 


0.98(8) 


-0.33(17) 


0.95 


<l>n 


1.84(6) 


1.20(08) 


-2.0(1) 


0.06 


1.59(5) 


-2.15(10) 


19.5 



TABLE II: Regge-plane fits combining both radial and angular-momentum trajectories (see main text for details). 



on the assumption that the probability of having a resonance 
with mass M and width F is of a Gaussian shape. 



Pi^/s) = Ce 



TM 



(Al) 



with C a normalization constant, whereas for the squared mass 
one has 



P{s)=C'e ^f2^^ 



(A2) 



The X^-fit then corresponds to applying the maximum- 
likelihood-method (MLM) and maximize with respect to a 



and Mq the function 



L(a,Mo;{r„,M„}) = nP(i„,r„,M„), 



(A3) 



n=l 



where s„ = an + Mq. This is the way the half-width rule is 
implemented in practice, i.e., by assuming short tails in the 
mass-distribution. On the other hand, from analyticity ar- 
guments the resonance profile function should be of a Breit- 
Wigner form, at least for sufficiently narrow resonances. Let 
us consider for definiteness the parameterization of a complex 
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resonance propagator at a given CM energy squared, s, 
Dis) = ^ 



s-M^-iTy/^' 



(A4) 



The ^/s in the denominator ensures that we have a pole on the 
second Riemann sheet (we neglect threshold effects). Like- 
wise, we also have a pure imaginary amplitude at the real res- 
onance value s = M^. The probability for such a mass distri- 
bution corresponds to the imaginary part ^, namely. 



PbwW=Z 



r^s 



(i-M2)2+r2i) 



(A5) 



where Z is a suitable normalization constant. Thus, we may 



apply the MLM to Eq. ( A3 1 for N resonances fulfilling the 



Regge formula and maximize with respect to a and Mq. 
As an specific example, to illustrate the difference between 



the Gaussian, Eq. ( A2 1, and the Breit-Wigner, Eq. ( A5 i, pro- 
files to the set of all 0++ scalars listed in the PDG (see also 
Fig. [8]) as discussed in Refs. [i33ii35J . where a joint formula 
for the trajectories. 



a 
-n- 

2 



(A6) 



was proposed. Maximizing P{a.ma;{Mn^Y„}) with respect 
to a and m<j and using Eqs. ( |A1| |, ( |A2[ ), and ( |A5[ ) yields the 
most likely values m^ = 0.545, 0.557, 0.562 GeV, and a = 
1.330, 1.336, 1.334 GeV^, respectively. The Gaussian cases 
correspond to the j^. analysis of Refs. 1331 l35l . Similarly to 
the ;|j 2. method, the errors can be determined by looking at the 

locus of the relative probability P{a,ma)/PmaK < {e^^^ '^), 
which for two variables yields A.x^ = 2.3 and 4.7 for the 68% 
and 95% confidence levels, respectively. We show the results 



in Fig. 18 where, as we can see, the resonance shape does not 
play a role. 

Of course, one may object to the previous confidence 
level analysis that for non-Gaussian probabilities mode (most- 
likely) and mean (average) are different. Indeed, the applica- 
tion of 

(A) = da dmaA{a,ma)LBw(a,ma;{M„,r„}) (A7) 

yields (a) = 1.34 GeV2 and (»%) = 0.53GeV for the mean 
values, whereas the mode is at m^ — 0.562 GeV and a = 

^ We are appealing to the Lehman representation for a resonance as obtained 
from a CM-energy dispersion relation of the scattering process, see, e.g.. 



1 .334 GeV2. Thus, for the Breit-Wigner case the mean and the 
mode are different and the errors are not defined by standard 
confidence level rules with (A) ± ^(A2) — (A)2, for which we 
get a = 1.34(8)GeV2, ma = 0.53(ll)GeV with a correlation 
r{a,ma) — —0.77. A way to sort this out is to define the equal 
probability contours, to integrate inside the inner region for a 
given confidence level. 



p{z} = dadmaP{a,ma)&{P{a,ma)—z), (A8) 



(N 



> 

o 

C3 1.3 



1.2 




m^ (GeV) 



FIG. 18: (color online) The 68% and 95% relative confidence-level 
contours of the radial Regge trajectories for all scalars in the PDG 
assuming several resonance profiles. 1) Gaussian at the level of the 
mass (dashed line), 2) Gaussian at the level of the squared-mass (dot- 
ted line) and 3) Breit-Wigner shape (black.solid). 



and to search for a z such that p{z) = 0.68. The resulting 
contour resembles strongly Fig. [18] reinforcing the conclu- 
sion that the shape of the resonance profile is irrelevant for the 
analyses of this work. 

Ref. 1401 for a discussion in the context of the ;r;r-scattering. 
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